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Abstract 

A first-principles-derived scheme, that incorporates ferroelectric and antiferrodistortive degrees of 
freedom, is developed to study finite-temperature properties of Pb(Zri_a;Ti3;)03 solid solutions near 
its morphotropic phase boundary. The use of this numerical technique (i) resolves controversies 
about the monoclinic ground-state for some Ti compositions, (ii) leads to the discovery of an 
overlooked phase, and (iii) yields three multiphase points, that are each associated with four 
phases. Additional neutron diffraction measurements strongly support some of these predictions. 
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The ferroelectric Pb(Zri_a.Tia;)03 (PZT) system is an example of perovskite solid solutions 
that are of high technological relevance because of their widespread use in piezoelectric 
transducers and actuators p|. Its phase diagram in a narrow compositional region centered 
around x=0.50, and denoted as the morphotropic phase boundary (MPB) [2], also makes 
PZT of large fundamental importance. For instance, recent measurements have discovered 
an unexpected ferroelectric (FE) monoclinic phase j3| - that leads to high electromechanical 
coefficients jj, ^ and acts as a structural bridge between the well-known FE tetragonal 
lases Slntl. Moreoever, another overlooked monoclinic phase has been 
Mft . wKich tKe „.„an.-oo™pet., R fe^oe.ect* disto.tion. 
antiferrodistortive (AFD) motions (that are associated with the oxygen octahedra rotation) 
coexist. However, many controversies still surround this latter monoclinic phase, such as 
its exact space group, the axis about which the oxygen octahedra rotate and even if it is a 



and rhombohedra 
recently observed 



minority phase rather than the long-range ground-state 

In view of these findings, one may also wonder if other surprises - such as other overlooked 
phase(s) or multiphase points (that are unusual thermodynamic features - can still be 
in store for PZT. Some features could indeed have been missed because of (i) the difficulties 
in growing and characterizing at various temperatures many samples with tiny compositional 
differences; and (ii) the current lack of theoretical tools that are able to accurately compute 
finite-temperature properties of perovskite alloys. For instance, the two sole first-principles- 
based finite-temperature schemes j^, 12| that yield a monoclinic phase in the phase diagram 
of disordered PZT have some shortcomings: the alloy effective Hamiltonian model of Ref. 

gives Curie temperatures that are ^ 60% larger than measurements and does not include 
AFD motions, while we are not aware of any study reporting the accuracy of the bond- 
valence model of Ref. for transition temperatures in PZT. 

The purpose of this Letter are three-folds: (1) to demonstrate that it is possible to develop 
a numerical scheme able to accurately compute finite-temperature properties of perovskite 
solid solutions and that incorporates both FE and AFD motions; (2) to apply such scheme 
to resolve the controversies discussed above, as well as, to reveal some striking overlooked 
effects in the MPB of PZT; and (3) to experimentally confirm some of these discoveries, via 
neutron diffractions. 

Our numerical scheme is based on the generalization of the first-principles-derived alloy 
effective Hamiltonian of Ref. to include AFD effects, in addition to FE degrees of freedom. 
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The total energy is thus written as a sumof two main terms, i?FE({ui}, {vh}-, {Vi}j Wj}) 
-Eafd-c({ui}, {'?/}, {o-j}, {cJi}), where ^^fe is the energy provided in Ref. p, while 

-Eafd-c is an additional term gathering AFD motions and their couplings in solid solutions. 
Uj is the local soft mode (which is directly proportional to the dipole) in unit cell i; {fju} and 
{rij} are the homogeneous and inhomogeneous strain tensors 1^, respectively; aj=+l or 
—1 corresponds to the presence of a Zr or Ti atom, respectively, at the B-lattice site j of the 
Pb(Zri_2,.Tia,)03 system. Finally, {wj} is a (B-centered) vector characterizing the direction 
and magnitude of the AFD motions in unit cell i. For instance, Ui = O.lz corresponds to a 
rotation of the oxygen octahedra by 0.1 radians about the z-axis. For Eafd-Cj we propose 
the following expression that contains five major terms: 

^AFD-c({u,}, {Vh}, {Vl}, Wj}, M) = + aAUjf + -iA^^X^^Xy + 'A,y^Xz + ^Xx^Xz)] 

i 

i,j a, 13 i a,f3 

i o,/3,7,(5 i j 

(1) 

where the sums over i run over all B-sites while the sums over j run over the six B 
nearest neighbors of the B-site i. rii{i) is the component, in Voigt notation, of the 
total strain (i.e., homogeneous and inhomogeneous parts) at the site i. a, /?, 7 and S 
denote Cartesian components, with the x-, y- and z-axes being along the pseudo-cubic 
[100], [010] and [001] directions, respectively. The first four terms of Eq. (1) characterize 
energertics of the hypothetical simple Pb< B >03 system resulting from the use of the 
virtual crystal approximation (VGA) [l4] (via averaging of the Zr and Ti potentials) to 
mimic Pb(Zro.5Tio.5)03, while the last term represents how the actual atomic distribution 
affects the AFD motions. The first four terms are the self-energy associated with AFD 
motions, the short-range interaction energy between AFD degrees of freedom, the interaction 
energy between strain and tUj, and the interaction energy between local- modes and AFD 
motions, respectively. The ka, ola and 7^ parameters, as well as the Kij^ap, Ci^ap and 
Dap'ys matrices (that greatly simplify due to symmetry) , are determined by performing local- 
density-approximation (LDA) calculations jl^ . Q] on 10 atom-supercells of Pb< B >03. 
On the other hand, the Aq and Ai parameters are evaluated from LDA calculations on 
supercells containing real Zr and Ti atoms in addition to VGA atoms. Note that the proposed 



expression for i?AFD-c differs from the one of Ref. [l^ for simple perovskites by its analytical 
expressions and the incorporation of alloying effects. 

The total energy of our effective Hamiltonian is used in Monte-Carlo (MC) simulations to 
compute finite-temperature properties of PZT systems. We use 12 x 12 x 12 supercells (8640 
atoms) and up to 4 million MC sweeps, as well as decrease the temperature in small steps, 
to get well-converged results. The {cxj} variables are chosen randomly to mimic disordered 
alloys. Outputs of the MC procedure identify the space group of the predicted phases, and 
are: the homogeneous strain tensor {r]H}, the < u > supercell average of the local mode 
vectors {u,}, and the < u >r vector defined as < u >r= jj ^. cijj(— l)"^W+"f (*)+"^(*), where 
the sum runs over the N sites i and where nx{i), ny{i) and nz{i) are integers locating the 
cell i (i.e., the B-site i is centered at [nx{i)^ + ny{i)y + nz{i)'z]a, where a is the lattice 
constant of PZT and where x, y and z are unit vectors along the Cartesian axes). A 
non- vanishing < u > indicates ferroelectricity while a non-zero < uj >r characterizes AFD 
motions associated with the R-point of the cubic first Brillouin zone. We also use the 
correlation functions of Refs to precisely determine transition temperatures. 

Figure 1(a) shows the x-, y- and z-Cartesian coordinates (< ui >, < U2 > and < >) 
of < u > in disordered Pb(Zro,52Tio.48)03 as a function of temperature, as predicted by 
the alloy effective Hamiltonian approach of Ref. that does not include AFD effects. 
Figures 1(b) and 1(c) display < u > and < lj >/j, respectively, in the same material, but 
when using the present method - with < uji >r, < uj2 >r and < 0^3 >r denoting the 
X-, y- and z-Cartesian coordinates of < a; >r. Figure 1(a) indicates that neglecting AFD 
motions and its couplings with ferroelectricity and strain leads to (1) a transition from a 
cubic paraelectric Pm3m phase (in which < u > is close to zero) to a tetragonal PAmm 
ferroelectric phase (in which < M3 > drastically increases while < ui > and < U2 > remain 
nearly null) around 825K ± 5K - which is higher than the experimental Curie temperature 
of ~ 663 ± 5 K 0] ; and (2) a tetragonal-to-monoclinic Cm transition, for which < ui > 
and < U2 > now increase and become equal to each other while still being smaller in 
magnitude than < M3 > , at 250K ± 13K. Comparing Figs l(b-c) with Fig 1(a) reveals 
that turning on AFD motions has dramatic effects on the phase transition sequence - as 
correctly guessed by Ref. [l8^. As a matter of fact, two new phases - with both FE and 
AFD components - are found at the lowest temperatures, namely monoclinic Cc between ~ 
163 K and ~ 138 K and rhombohedral i?3c below 138 K ± 6 K. Furthermore, the theoretical 
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Curie temperature is decreased down to 413 K ± 13 K, which now makes it smaller than 



the measurement 17[. This discrepancy is due to the fact that the effective Hamiltonian 
parameters were determined at the (underestimating) theoretical LDA lattice constant. As 
a matter of fact, we numerically found that the pressure leading to our theoretical volume 
matching the experimental one at low temperature is —4.68 GPa, and that the use of 
this negative pressure within the present numerical scheme results in the < u > and < uo >fi 
displayed in Figs 1(d) and (e), respectively. The resulting theoretical Curie temperature is 
now ~ 650 ± 13 K, which is in excellent agreement with the measured value of ~ 663 ± 
5 K. Moreover, the symmetry of the purely FE phases coincides with those measured at 
this specific composition |6], that is tetragonal PAmrn and monoclinic Cm when decreasing 
the temperature. Our simulations also predict the existence of a monoclinic state below 
138 K ± 13 K in which ferroelectricity cohabits with AFD motions, as consistent with the 
electron-diffraction reflections observed in Ref. below 150 K. Furthermore, we confirm 
the findings of Refs 0, Q, 0] that the correct space group for the low-temperature phase of 
Pb(Zro.52Tio.48)03 is Cc rather than Pc, as initially proposed in Ref. [9]. Moreover, this 
Cc phase is indeed the (long-range-ordered) ground-state of Pb(Zro.52Tio.48)03 rather than 
a minority phase coexisting with the purely FE Cm phase for the lowest temperatures, as 
suggested in Ref. [^. Finally, the oxygen octahedra in this Cc phase do neither TotaXe about 
;he [001] direction (as thought in Ref. 0]) nor about the [111] direction (as proposed in Ref. 
^) but rather about an axis that is in-between the [001] and [111] directions (which is in 
agreement with Ref. Q]), since Fig 1(e) shows that, at the smallest temperatures, < c^a >r 
is larger in magnitude than < cji >/? and < uji >r= < uj2 >b- Interestingly, this axis of 
rotation does not precisely coincide with the direction of the polarization - since the "FE" 
ratio differs from the "AFD" ratio ^"^^^^ - which contrasts with the assumption of 



Ref. 

We now turn our attention to Fig. 2 that displays the phase diagram of Pb(Zri_a;Ti2,)03 
near the MPB, as predicted by the present scheme with a —4.68 GPa pressure. Seven phases 
exist within a narrow compositional range. Five of such phases are well-known, namely para- 
electric Pm3m, FE rhombohedral R3m, FE tetragonal P4mm, FE monoclinic Cm, and the 
rhombohedral R3c state in which ferroelectricity coexists with AFD motions. The sixth 
phase is the previously "controversial" Cc phase discussed above, and that is the ground- 
state at intermediate Ti compositions. (Note that the monoclinic phases are predicted to 
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exist for the narrow Ti compositional range of 46.3 %-51.5%, which agrees remarkably well 
with the experimental range of ~ 46%-52% 19 1.) The last phase occurs at small tem- 
perature and for the largest displayed Ti compositions, and has never been suggested to be 
a possible ground-state of PZT (to the best of our knowledge). This phase has a /4cm space 
group, also involves a coexistence of ferroelectricity and rotation of oxygen octahedra, but 
is associated with the tetragonal symmetry (unlike R3c and Cc). As a result, an increase in 
Ti concentration from ~ 47% to 52% results not only in a spontaneous polarization contin- 
uously rotating from [111] to [001] but also to the oxygen octahedra varying their rotation 
axis from [111] to [001] - via the [uuv] directions in the bridging monoclinic structures Cm 
and Cc phases. [We also checked that no AFD state associated with an off-center k-point 
different from R exists in the MPB of PZT.] Other features of Fig. 2 going against common 



beliefs p, |8| are the non-verticality of the temperature-composition boundary lines separat- 
ing rhombohedral from monoclinic structures. As a result, five different phases emerge for 
some narrow compositions. For instance, for x ~ 47.2%, our scheme predicts the Pm3m, 
PAmm, Cm, R3m and R3c phases when decreasing the temperature. Similarly, a Ti con- 
centration around 47.5% yields a Pm3m - PAmm - Cm - Cc - R3c sequence 20]. Finally, 
Fig. 2 further reveals the existence of three temperature-compositional points, for which 
four phases are very close to each other! (Multiphase points are unusual thermodynamic 
features Q].) The first multiphase point gathers the paraelectric Pm3m state and the FE 
i?3m. Cm and PAmm phases, and occurs for T ~ 610 K and x ~ 46.3 % (Note that 

a quadruple point in which the paraelectric state coexists with three FE phases was also 
predicted to exist in the phase diagrams of epitaxial PZT thin films j3].) Two purely FE 
phases meet with two phases exhibiting a coexistence of FE and AFD motions for the other 
two multiphase points. More precisely, i?3m. Cm, R3c and Cc all meet at T ~ 150 K and 
X ^ 47.4%, while Cm, PAmm, Cc and /4cm overlap at T ^ 100 K and x ^ 0.495. 

We have also undertaken a Rietveld analysis of neutron data on a Pb(Zro.49Tio.5i)03 
sample, collected at lOK at the ILL thermal source on the powder diffractometer D2B 
(A=1.6 A). Analysis of the data collection from 2G=20 to 140 degrees with a 2B step of 

231 ] ■ using pseudo-Voigt peak-shape 



0.05 degrees was carried out with the XND software 
function including asymmetric broadening and linear interpolation for the background. The 
fluctuation in the sample's composition was estimated to be ~ 0.003 based on x-ray peaks 
at 300K and dielectric measurement. We observed superlattice reflections in the diffraction 
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pattern (not shown here), and thus considered different space groups that are associated with 
doubled unit cells - namely, rhombohedral i?3c, monoclinic Pc, Cc and /c, and tetragonal 
lAcm. Each single-phase model leads to a considerable mismatch between the observed and 
calculated profiles for the pseudocubic reflections, and thus to poor agreement factors. A 
considerable improvement is achieved by considering two-phase coexistence models between 
the above space groups and also by including the monoclinic Cm and cubic Pm3m phases. 
The best agreement factors and matching between the observed and calculated profiles is 
obtained with the lAcm+Cc coexistence model (that results in R^p, and G.o.F. factors 
being equal to 6.69%, 5.09% and 1.72, respectively). Our structural refinement thus confirms 
the existence of the "controversial" Cc state (for which our structural parameters are very 
close to the ones of Ref. j^]) and of the presently-discovered /4cm structure (for which the 
a, b and c parameters are equal to 5.700(2), 5.700(2) and 8.329(6) A, respectively). Note 
that the proportion of /4cm and Cc phases is approximately 7:3. The coexistence of these 
two phases may be due to compositional fluctuation in the sample, which will make this 
coexistence consistent with the theoretical phase diagram of Fig.2. It rnay also take its 
origin from the presence of small grains in the grown ceramic, since Ref. j2^ predicts that 
the coexistence of different phases may be energetically favorable (over a single phase) in 
a ceramic of fixed composition when the grain size is small. We also hope that our other 
predictions, including the ones about multiphase points, will be experimentally confirmed 
soon. 
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thank B. Noheda for useful discussions. 



[1] K. Uchino, (Kluwer Academic Publishers, Boston) (1996). 

[2] M.E. Lines and A.M. Glass, Principles and Applications of Ferroelectrics and Related Materials 

(Clarendon Press, Oxford, 1977). 
[3] B. Noheda, D.E. Cox, G. Shirane, J. A. Gonzalo, L.E. Cross, and S-E. Park, Appl. Phys. Lett. 

74, 2059 (1999). 

[4] B. Notieda, D.E. Cox, G. Shirane, S.E. Park, L.E. Cross, and Z. Zhong, Phys. Rev. Lett., 86, 



7 



3891 (2001). 

[5] L. Bellaiche, A. Garcia and D. Vanderbilt, Phys. Rev. Lett. 84, 5427 (2000); Ferroelectrics 

266, 41-56 (2002). 
[6] B. Noheda, L. Wu, and Y. Zhu, Phys. Rev. B 66, 060103(R) (2002). 

[7] D.M. Hatch, H.T. Stokes, R. Ranjan, Ragini, S.K. Mishra, D. Pandey and B.J. Kennedy, 

Phys. Rev. B 65, 212101 (2002). 
[8] D.J. Woodward, J. Knudsen and I.M. Reaney, Phys. Rev. B 72, 104110 (2005). 
[9] R. Ranjan, Ragini, S.K. Mishra, D. Pandey and B.J. Kennedy, Phys. Rev. B 65, 060102(R) 

(2002). 

[10] D. Vanderbilt and W. Zhong, Ferroelectrics 206, 181 (1998). 

[11] M. L. Plumer, A. Caill and K. Hood Phys. Rev. B 39, 4489 (1989). 

[12] I. Grinberg, V.R. Cooper and A.M. Rappe Nature 419, 909 (2002). 

[13] W. Zhong, D. Vanderbilt and K.M. Rabe, Phys. Rev. Lett. 73, 1861 (1994); Phys. Rev. B 52, 
6301 (1995). 

[14] L. Bellaiche and D. Vanderbilt, Phys. Rev. B 61, 7877 (2000). 

[15] P. Hohenberg and W. Kohn, Phys. Rev. 136, B864 (1964); W. Kohn and L.J. Sham, ihid. 
140, A1133 (1965). 

[16] A. Garcia and D. Vanderbilt, Appl. Phys. Lett. 72, 2981 (1998); First-Principles Calculations 
for Ferroelectrics: Fifth Williamsburg Workshop, R.E. Cohen, ed. (AIP, Woodbury, New York, 
1998), p. 53. 

[17] B. Noheda, J. A. Gonzalo, L.E. Cross, R. Guo, S.-E. Park, D.E. Cox and G. Shirane, Phys. 

Rev. B 61, 8687 (2000). 
[18] M. Fornari and D.J. Singh, Phys. Rev. B 63, 092101 (2001). 

[19] K. C. Lima, A. G. Filho, A. P. Ayala, J. Mendes Filho, P. T. Freire, F. E. Melo, E. B. Arajo, 
and J. A. Eiras, Phys. Rev. B 63, 184105 (2001). 

[20] Ref. ^ indicates that the Ti composition, for which the PAmm-to- Cm transition occurs at 
300K is ~ 48% while our scheme predicts that this happens for x ~ 47.3%. This slight dis- 
agreement may be due to the difficulty for the simulations to accurately reproduce the strong 
compositional dependency of the PAmm-io-Cm transition temperature (see Fig. 2), or that 
real samples have compositional inhomogeneities. 

[21] The simulations in Fig. 2 provide Curie temperatures that are in very good agreement with 

8 



measurement for x ~ 48%, but that begin to underestimate (respectively, overestimate) ex- 
)erimental data by more than 50 K for x smaher than 46.5% (respectively, larger than 50.5%) 
^. As a result, the first multiphase point may be experimentally found at a temperature 
slightly larger than 610 K. 

[22] N.A. Pertsev, V.G. Kukhar, H. Kohlstedt and R. Waser, Phys. Rev. B 67, 054107 (2003). 

[23] J.F. Berar, lUCr. Sat. Meeting on Powder Diffractometry, Toulouse (1990). 

[24] R. Ranjan, A.K. Singh, Ragini, and D. Pandey, Phys. Rev. B 71, 092101 (2005). 

[25] A.G. Zembilgotov, N.A. Pertsev and R. Waser, J. Appl. Phys. 97, 114315 (2005). 



9 



FIG. 1: Supercell average < u > of the local mode vectors (Panels (a), (b) and (d)) and AFD- 
related <uj >r quantity (Panels (c) and (e)), as a function of temperature in Pb(Zro.52Tio.48)03. 
Panel (a) corresponds to predictions from the use of the alloy effective Hamiltonian method of 
Refs. 5| - that neglects AFD effects. Panels (b) and (c) (respectively, Panels (d) and (e)) shows 
the results of the presently proposed method with no applied pressure (respectively, when applying 
a pressure of — 4.68GPa). 

FIG. 2: Phase diagram of Pb(Zri„a;Ti^)03 near its MPB, as predicted by the present scheme with 
an applied pressure of — 4.68GPa. Symbols display the direct results of our simulations, while lines 
are guide for the eyes. Indices 1, 2 and 3 indicate the multiphase points. The uncertainty on the 
transition temperatures is typically around 13 K, except close to the multiphase points 2 and 3 for 
which this uncertainty is around 3K. 
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